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ABSTRACT 



We show that the use of Doppler shifts of Zeeman sensitive spectral lines to 
observe waves in sunspots is subject to measurement specific phase shifts aris- 
ing from, (i) altered height range of spectral line formation and the propagating 
character of p mode waves in penumbrae, and (ii) Zeeman broadening and split- 
ting. We also show that these phase shifts depend on wave frequencies, strengths 
and line of sight inclination of magnetic field, and the polarization state used for 
Doppler measurements. We discuss how these phase shifts could contribute to 
local helioseismic measurements of 'surface effects' in sunspot seismology. 

Subject headings: Sun: helioseismology — Sun: magnetic fields — Sun: oscilla- 
tions — radiative transfer — sunspots 



1. Introduction 



The solar oscillations are measured with better signal-to-noise from Doppler shifts of a 
spectral line than intensity fluctuations. Existence of reliable algorithms to estimate Doppler 
shifts from lirie inte nsities (filtergrams) scanned over a few locations across a line (MDI; 



Scherrer et al.l (Il995l )). or a fast Fourier tachometric scanning of a line (GONG: lHarvey et al. 



(I1994J )). makes possible deriving velocity images of solar surface frequent enough to study 
the solar oscillations. Above methods of estimating line of sight velocities from Doppler 
shifts assume rigid displacements of a spectral line in wavelength in response to material 
motions. However, gradients along line of sight in solar velocity and other physical quantities 
introduce asymmetries in line profile leading to Doppler cross-talk of signals from different 
layers in the atmosphere. In quiet Sun, the five minute band p mode oscillation signals are 
evanescent at observation heights and hence oscillation phases are only minimally affected 
by height gradients. In active regions, not only the vastly different thermal and magnetic 
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structure complicate Doppler shift measurements through changes in hne profile shapes but 
also altered character of waves in line forming layers. 

In this Letter, using a temporal sequence of spectropolarimetric observations of a sunspot 
we demonstrate: (i) propagating character of p modes in the penumbra and extended height 
range of spectral line formation cause phases of waves to depend strongly on intensity level 
within a line from where Doppler velocities are measured, (ii) sampling Zeeman split portions 
near the line core causes filtergraph (MDI/SOHO) or bisector methods for Doppler velocity 
estimation to introduce spurious phase shifts in the waves measured. Wave propagation 
signatures, in the firs t case above, which have recently been studied in chromospheric lines 
(IFinsterle et al.ll2004l ). are also shown to differ in the Doppler shifts of Stokes I and polarized 
components (circular and linear). We discuss how the above findings are related to recen t 
local helioseismic measurements of 'surface effects' (jSchunker et al.ll2005l ; IZhao et al.ll2006l ). 



2. Data and Analysis Methods 



The observations were carried out using the Advanced Stokes Polarimeter (ASP) at the 
Dunn Solar Telescope of the National Solar Observatory at Sacramento Peak, Sunspot, New 
Mexico. A medium sized sunspot (NOAA AR0750, diameter ^ 16 Mm) located close to disk 
center was observed on April 12, 2005, in mid-photospheric spectral line Ni I (A=6768 A, 
used by MDI/SOHO and GONG). We have made a temporal sequence of Stokes vectors of 
this line by scanning the sunspot over consecutive (in x-direction) 14 slit positions oriented 
N-S (y-direction) in one minute and repeating it for about 100 minutes (with a solar image 
tracker fixing the field of view). The N-S extent of slit covers about 64 Mm with a resolution 
of 0.375" with the spot at the center and, with a slit width of 0.6", the 14 slit positions 
(x-direction) cover mainly umbral region. 

We do Milne-Eddington (M-E) inversions of Stokes vectors (/, Q, U, V^s) to obtain line 
of sight velocities v{x,y,t), total magnetic field strengths B{x,y,t), line of sight inclination 
jix,y,t) and azimuth 'iljix ,y,t) of magnetic field among other standard inverted quantities 
(ISkumanich fc Lited 119871 ). The M-E inversion procedure is weighted towards polarized 
component within magnetic field a nd provides velocity estirnates from the line co re position 
which is used as a free parameter (ISkumanich fc Lited 119871 : IJefferies et al.lll989l ). Since we 
want to examine how wave phases change with height within line formation layers, how they 
are seen in polarized and unpolarized light, and how they contribute to local helioseismology 
measurements of 'surface magnetic effects', we also obtain line of sight velocities from the 
fol lowing methods: (i) Doppler shifts of bisector points at diffe rent intensity levels within 



/ (IKeil fc YackovichI Il98ll : ICavallini et al.l Il985l : iRimmeld Il995l ) , and those within circular 



- 3 - 



(average of I + V and I — V velocities, hereafter denoted CP) and linear (I + Q, denoted LPl) 
polarized components (ISankarasubramanian fc Rimmeld l2002l : Idel Toro Iniesta II2003I ). and 
(ii) the two Doppler filtergram measur ements in circular ( CP) and linear (LPl) polarized 
components employed by MDl/SOHO (IScherrer et al.lll995l ). 



Phase shifts between any two velocity signals, fi(t) and V2{t), are calculated according 



to, 



tan 



1 / Im[vi{ uj)v;{uj)] 
\Re[vi{. 



where Vi{uj) and V2{u)) are the Fourier transforms of the two time series. In the above 
convention, positive values for 0i_2 at positive u mean vi(t) is advanced in phase with respect 
to V2{t). Our aim here is to study how (p between velocities from different bisector levels as 
well as that between different measurements vary over the sunspot as a function of position 
and u, i.e. (j){x,y, z,u>), due to propagating waves and altered line profile shapes (Zeeman 
splitting). The vertical co-ordinate z represents height within line forming layers and we 
use intensity levels (i.e., bisector points) within the line to scan it. Based on the plane- 
parallel VAL-C model of quiet solar atmosphere, the Ni I line formation is reported to span 
a height range of 18 km (continuura) - 288 km (line c ore) {z=0 being referenced to continuum 
optical depth Tc=1) (INorton et al.ll2006l : lBrulslll993l ). We do not attempt full Stokes profiles 
inversions to determine z dependences of this is not only computationaly expensive (as we 
have too many profiles from the x, y, t scans of Stokes vector) but also lead to noisier signals 
than those from direct bisector analyses. In quiet-Sun region the main p mode band acoustic 
waves (a; < oOac ~5.3 mHz, the acoustic cut-off frequency) are largely reflected down into the 
solar interior, and hence (p{x,y,z,LL!) ^ constant = 0; at a; > Uac, however, because of wave 
propagation, 0(x, y, z, uj) ~ ^(z, u) ^ 0. Within the sunspot, we study x and y dependences 
of using temporal averages of M-E inverted B[x^ ?y, t) and 7(x, y, t) shown in Figure 1. 



3. Surface magnetism effects due to propagating waves and Zeeman splitting 

In order to analyse phase shifts due to the propagating nature of waves, we calculate 
bisector velocities from 9 equally spaced intensity levels within / and CP profiles, Ji=0.1Jc, 
J2=0.2/c,..., /9=0.9/c, where Ic is the continuum level. Using the time series of velocities from 
different bisector levels, we calculate phase shifts between them using Equation [H When 
the line is Zeeman split, formation heights of bisector levels of / profiles near the core are 
significantly altered and any asymmetries in depths of Zeeman components further corrupt 
the velocity signals. Bisector points of individual CP components are largely free of the 
above uncertainties. In view of the above, we use M-E inverted velocities as representative 
of the highest layer within line formation region, and we show that this is indeed the case 
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by comparing them with CP bisector velocities from near the core for the Ni I hne. 

The phase shifts (f){x,y, z,uj) exhibit random fluctuations over space as well as in a; due 
to the inherent stochastic nature of the solar velocity field and due to noise from various 
sources. Fluctuations over cu are reduced by taking median values of (j) over the p mode 
band (2-5 mHz) or over bands of 1 mHz width centered at 3, 4 and 5 mHz (to study any 
frequency dependence). Signals over space are studied using the temporal averages B{x,y) 
and j{x,y) (Figure 1): we average the phase shifts over all pixels that fall within a 25 G 
window of B or 2° of 7 and assign the average to all those pixels (spatially) and to the 
central values of B and 7. The resulting spatial maps of (j) arc essentially one dimensional 
functions of 5 or 7 with any azimuthal variations due to such variations in S or 7 averaged 
out. In Figure 2, we plot (j) between mid-level {I5 — 0.5/c) bisector velocities and those from 
other levels: panels a and b show median (f)^f{i— 1,3,7,9) (superscripts identify the profiles 
used and subscripts the pair of levels that the phase shifts come from) over the full p mode 
band as a function of B and 7, respectively, panel c shows ^^f and (p^g as a function of 
B for different frequency bands, and panel d compares (j)'"^ and (j)^ between levels (5,1) and 
(5,9) as a function of B for the full p mode band. For clarity, we have plotted error bars 
only at representative values of B or 7. The error bars correspond to standard deviations 
of (j) over the pixels that fall within the bin sizes of B or 7. The spatial maps that are 
obtained after binning over 25 G windows in line of sight magnetic field B^ — Bcos{'-f) are 
shown in Figure 3: in panels a and b, respectively, are ^5 ^ and </'^f , where i =M-E,1,2 
and 9, in panel c are 0^'^^ that are between same level bisector velocities from GP and / 
profiles, and in panel d are (p^ g and (f){ g in the 2.5 - 3.5 mHz and 3.5 - 4.5 mHz frequency 
bands. Ghanges in against B or B^ and the bisector intensity levels, in Figures 2 and 3, 
clearly show that these are due to propagating waves, even within the main p mode band, 
and that propagation is mainly confined to penumbral region. In quiet Sun pixels {B < 
200 G), due to wave propagation into higher layers at a; > ujac, though much smaller 
than that seen in the penumbra, is noticeable in panels c and d of Figure 2. Both Stokes / 
and GP profiles' velocities show propagating waves but the signals are larger in GP bisector 
velocities indicating that propagation is mostly in magnetized channels (panel d of Figure 
2 and panels a and b of Figure 3). Phase shifts between same level bisector velocities from 
I and GP profiles, shown in panel c of Figure 3, further substantiates the above inference. 
The 0^'*-^^ from wing level (7 and 9, the right two panels in panel c of Figure 3) bisector 
velocities are very small, suggesting that differences in formation heights of wings of / and 
GP profiles is small. 

The above results show that if we have to minimize contributions from propagating 
waves we should restrict Doppler measurements to wings. In panel d of Figure 2 and panel 
a of Figure 3, it is seen that the Stokes I velocities introduce spurious advanced phases 
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within the umbral and umbral - penumbral boundary region, when bisector velocities from 
near the core are involved: the fallen or collapsed 05 ^ plotted as red triangles in panel d of 
Figure 2, and the scrambled phase signals seen in the middle two maps of panel a in Figure 
3. We attribute the above to changes in / profile shapes due to Zeeman splitting, as these 
signals are not present or weak in CP measurements. We confirm this using the M-E inverted 
velocities, which model the Zeeman splitting and retrieve velocities; the left most map in 
panel a of Figure 3 show that 4>i^ME free of the above discussed artifacts, and also matches 
well with map of (P'^me panel b of Figure 3. The similar values and structure of 



and (1)^2 to that of (p'^ME further clarifies that wave propagation is confined to magnetized 



kCP 



'5,ME and (pl%jE m 



medium. Interestingly, the M-E velocities show negative changes for 
the outer penumbra (panels a, b of Figure 3), which we identify as a possible signature of 
downward propagating waves along field lines that return back to the photosphere. 



The MDI measurements (jScherrer et al.lll995l ) use CP and LPl profile s: most full-disk 
observations are based on LPl and almost all he l ioseismic holography resu lts (jSraun fc Lindsey 
I999I : iBraun et al.l 120041 : ISchunker et aP boosi : llindsev fc BraunI boOSal lbl) on surface mag- 
netism effects are based on these data; on t he other hand, most tin i e-distance heliose ismology 
studies of sub-surface struc ture and flows (Kosovichev et al. 2000; Zhao et aD 2001 ) and the 
surface magnetism effects ( Zhao et al. 20061 ) have used the MDI CP hi-res data. Here, we 
have derived Doppler velocities from both of these MDI measurements by executing the MDI 
algorithms (which the onboard processor of MDI uses). In quiet Sun, i.e. when there are no 
strong magnetic fields to modify the line profile through Zeeman effect, both of these MDI 
measurements reduce to deriving Doppler velocities from unpolarized Stokes / and hence 
are the same. However, as demonstrated in Figures 2 and 3, propagating waves introduce 
phase shifts, and to estimate how much of these are present in MDI measurements we have 
calculated phase shifts between MDI CP and LPl velocities and wing level (level 7) bisector 



velocities from /, denoted respectively as 



^Ifoi and 



lLPI 
'MDI 



. The results are in Figure 4: panel 



a shows and panel b shows 0^^/ over bands of 1 mHz width centered at frequencies of 

3, 4 and 5 mHz and also the full p mode band (2-5 mHz). We found from bisector velocities 
from near the core that Zeeman splitting causes an advancement of phases within umbra 
(panel a and d in Figure 3); since LPl profile within umbra is the same as / (because 7 is 
close to zero), similar advanced phases are expected from MDI LPl measurements. This is 
clear in panel a of Figure 4. In summary, for the medium size spot observed, surface effects 
due to wave propagation in the penumbra cause a phase lag in the range of 5 - 10 seconds 
in both MDI LPl and CP measurements and the Zeeman split profiles cause positive phase 
shifts of the order of 8 seconds in umbra in MDI LPl measurements. 
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Discussions and conclusion 



Recent stud ies in sunspot seismology ( Schunker et al. 2005 : Lindsey &: Braun 2005al lbl: 
Zhao et al.l l2006l ) have recognised the role of 'surface magnetic effects' that contribute to 
seismically measured phases or travel times of acoustic waves. We have shown that phases 
of acoustic waves observed within sunspots using Doppler shifts of a spectral line have 
signatures of physical changes that waves undergo within line forming layers as well as of 
systematics in Doppler measurements induced by Zeeman split profiles. We have also shown 
that wave propagation effects and systematics from split portions of the line profiles are 
minimal when Doppler velocity measurements are restricted to wings of a spectral line. 
Since a seismologically correct determination of phase shifts require separating the above 
surface effects in the measurements, our latter results above show that such a correction 
can be carried out in the observation procedure itself: within sunspots, the Doppler shift 
measurements could be restricted to wings of spectral lines. In the case of MDl/SOHO, 
the onboard algorithm combines all the filtergrams to determine Doppler velocities. Since 
the filter positions are fixed with respect to the central wavelength, when the line is wider 
due to Zeeman splitting, there is a systematic sampling of split core region of the profiles 
and hence higher atmospheric layers. This causes phase lags in waves from penumbral 
region and spurious advanced phases within umbra when LPl profile is used (most full-disk 
measurements). Correcting the MDI measurements for such surface effects, hence, requires 
maps such as the ones derived here (Figure 4). A detailed 'calibration' of phase shifts due 
to surface effects against B or j require similar studies involving spots of various sizes. 

We have not assessed how the surface effects studied here would depend on positions 
across the solar disk. However, the present finding that propagating waves in penumbra cause 
a phase shift depending on height of observation points to similar effects when there are line 
of sight changes in optical depth. Such changes are expected for sunspots located away from 
disk center, due to the Wilson depression, and this combined with the field line alligned wave 
propagation can introduce opposite changes, in measured phases, in the limb side (deeper 
level of observation) and center side (higher level of observation) penumbrae. We suggest that 
such effects on the phase shifts are a lik e ly con tribu tor to the azimuth ally varying 'inclined 
magnetic field effect' of ISchunker et al.l (120051 ) and IZhao et al.l (120061 ). The magnitude of 
surface effect signal that we have, for the medium sized spot, falls in the range of 5 - 15 
seconds (Figures 2, 3 and 4) and it points to larger phase shifts for larger spots. We note 
that, for fairly large sunspots, the m agnitude of travel time perturbations measured in time- 



distance helioseismic measurements (IDuvall et al.l Il996l : iKosovichev et al. 



200G 



Zhao et al. 



20011 ). as well as the vantage dependent control correlation phases (ISchunker et al.ll2005l ). fall 
in the range of 30 - 60 seconds. Our results on surface effect signal due to Zeeman splitting 
in umbral region (Figures 3 and 4) call for a careful study of active regions in one or more 
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magnetically insensitive photospheric lines. 
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Fig. 1. — Maps of continuum intensity Ic and temporal averages of M-E inverted B and 7 
over the sunspot from Ni I observations. 
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Fig. 2. — Phase shifts between velocities from different bisector levels of CP and / profiles 
of Ni I line: panels a and h show (fP^ against B and 7, respectively, for the full p mode 
band (2-5 mHz), panel c shows ^^f and 0^1" over different frequency bands, and panel d 
compares (j/"^ and (f)^ from the bisector levels (5,1) and (5,9) for the full p mode band; see 
text for further details 
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Fig. 3. — Maps of phase shifts of waves observed using Ni I hne; see text for details 
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Fig. 4. — Maps of phase shifts due to 'surface magnetic effects' in MDI hnear (LPl, full disk) 
and circular polarization (CP, hi-res) measurements. 



